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Abstract: We argue that the solar g-modes are unlikely to have caused the discrete peaks 
in the power spectrum of the solar wind flux observed by Thomson et al. (1995). The 
lower limit to the energy of individual g-modes, using the amplitudes given by Thomson et 
al., is estimated to be at least 10^^ erg for low order g-modes; the resulting surface velocity 
amplitude is at least 50 cm s~"^, larger than the observational upper limit (5 cm s~"^). 

We suggest that the most likely source for the excitation of solar g-modes is turbulent 
stresses in the convection zone. The surface velocity amplitude of low degree and low 
order g-modes resulting from this process is estimated to be of order 10~^ cm . This 
amplitude is interestingly close to the detection threshold of the SOHO satellite. The long 
lifetime of g-modes (~ 10^ years for low order modes) should be helpful in detecting these 
small amplitude pulsations. 
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1. Introduction 

Gravity mode oscillations of the sun are primarily confined to its radiative interior and 
their observation would thus provide a wealth of information about the energy generating 
region which is poorly probed by the p-modes. In the past 20 years a number of different 
groups have claimed to detect g-modes in the sun (e.g. Brookes et al. 1976; Brown et al. 
1978; Delache & Scherrer 1983; Scherrer et al. 1979; Severny et al. 1976; for a detailed 
review of the observations please see the article by Palle, 1991, and references therein), but 
thusfar there is no consensus that g-modes have in fact been observed. Recently Thomson, 
MacLennan & Lanzerotti (1995), hereafter referred to as TML, have reported detection 
of g-modes in the flux of high energy particles in the solar wind. In the absence of a 
detailed model of how pulsation affects the solar wind flux, it is diflicult to determine the 
velocity amplitude of the pulsation in the photosphere from the amplitude of the solar 
wind flux variation. We, however, use general physical arguments to make an estimate of 
the photospheric velocity amplitude of g-modes from the wind data (section 2) which can 
be compared with the results of optical searches. In section 3 we provide estimates of the 
surface velocity amplitude, arising from excitation by turbulent convection, of g-modes of 
several different degree [£) and order (n) for comparison with this and future observations. 
The main results are summarized in section 4. 



2. Observation of g-modes in the solar wind? 

The dispersion relation for solar g-modes in the limit of n ^ ^ is (e.g. Christensen- 
Dalsgaard & Berthomieu, 1991) 

P,, = , ^° (2n + e-5) 



where 

Nb is the Brunt- Vaisala frequency, ri is the upper turning point of the mode, and 5 5/6. 
Thus the order n of a g-mode of frequency 5 /iHz (the low frequency end of the observed 
peaks attributed to g-modes by TML), and degree ^=1,2, and 3 is about 130, 230, and 
320, respectively. Thomson et al. flnd several groups of very closely spaced frequencies 
in their data, for instance 106.768, 111.938 and 114.042 /iHz (see Table 2 of TML). These 
can not be modes of the same ^, if indeed they are g-modes, since the observed frequency 
spacing is considerably smaller than the frequency spacing expected for modes of flxed ^, 
but larger than that due to rotational splitting (high degree g-modes have smaller frequency 
spacings, but they are unlikely to be observed because their surface velocity amplitudes are 
very small due to large attenuation in the convection zone). Thus, if g-mode amplitudes 
in the solar wind have a frequency dependence similar to that of table 2 of TML, then one 
expects there to be about 2000, £=1 & 2, g-modes with frequencies in the observed range. 
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which should have amphtudes of order those detected by TML. (This is a conservative 
estimate since TML attribute some of their peaks to modes of ^ > 2.) Thomson et ah, 
however, detect only about 40 such modes, which is a surprisingly small fraction of the 
total number of observable modes. Although there may be unstated reasons why other 
frequencies were not observed, it seems unlikely from an a priori point of view that such a 
small subset of frequencies should be detected. 

The mean mass flux in the solar wind is ~ 1.3x10"^^ g (or 2xl0~^'^ Mq per year), 
and the mean wind speed is 400 km . Thus the total kinetic energy flux in the wind is 
about 10^^ erg s~"^; most of this energy is carried by the thermal particles, while the high 
energy particles observed by TML carry several orders of magnitude less energy. TML 
flnd the fractional variation in the high energy particle flux in the solar wind associated 
with g-mode pulsation to be about 2% for low frequencies (~ 10/i Hz) and 0.3% for high 
frequencies (~ 100/i Hz). Since the fractional variation in the magnetic fleld is also ~2% 
(Thomson, personal communication), the flux of thermal particles in the wind is expected 
to be modulated by an amplitude similar to that of the high energy particles. Thus, 
provided that a signiflcant fraction of the flux variation of the thermal particles is due to a 
change in the wind speed, TML's observations imply that the g-modes lose energy to the 
solar wind at a rate of about 10^"^ erg for low £ modes with n > 10 and 10^^ erg for 
modes with n of order a few (the mean rate of energy loss is proportional to the square of 
the variation in the particle flux). An independent estimate of the rate of g-mode energy 
loss can be made using the observed variation of the magnetic fleld strength associated with 
a mode in the solar wind. These measurements are made using near-Earth satellites and 
are also reported by TML. Magnetic fleld variations are expected to propagate as Alfven 
waves, which carry energy away from the sun at a rate of ~ c/^(^i?)^ V^/2, where d is the 
distance of the spacecraft from the sun, SB ~ 3xl0~^ gauss is the magnetic fleld variation 
associated with a mode (Thomson, personal communication), and Va ~ 50 km is the 
Alfven speed in the solar wind. Thus, we infer that a g-mode loses approximately 5x10^"^ 
erg to Alfven waves in the wind, which is consistent with the above estimate for low 
n modes. 

If we assume that the damping of the g-modes due to their modulation of the solar 
wind is not the dominant damping mechanism for these modes, then the lower limit to 
the g-mode energy is estimated to be about 10^^ erg for low £ modes with n of order a 
few and 10^^ erg for low £ modes with n > 10. These estimates were obtained by making 
use of the dissipation time (turbulent plus radiative) of g-modes, which is of order 10^ 
years for low order modes and decreases with n as (n is the mode order). The resultant 
surface velocity amplitude is at least 2 cm for low £ modes with n > 10 and 50 cm 
for modes with n of order a few. The latter amplitudes are larger than the observational 
upper limit of 5 cm (Kuhn, Libbrecht and Dicke (1986) flnd a limit, for ^ > 2, of about 
15 cm at 100 /xHz and about 2 cm at 10 /iHz; Garcia, Palle and Cortes (1988) 
flnd a limit of about 4 cm for low £, including £=1; Frohlich (1990) also flnds a limit 
of about 4 cm s~"^; for a review of the results of various optical searches for g-modes, see 
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Palle, 1991, and references therein). 

Thomson et al. also find p-modes in the solar wind electron flux data. The fractional 
variation of the electron flux associated with a 3 mHz p-mode is about 10""^ (Thomson, 
personal communication). Using the argument of the preceding paragraphs, with a mode 
lifetime of 10 days, we estimate the lower limit to the p-mode surface velocity amplitude 
to be about 1 cm s~"^; the observed amplitude is about 10 cm . 

Particle flux perturbations imposed at the solar surface are dispersed out over a dis- 
tance of a few tens of solar radii due to the distribution of particle speeds (this is damping 
due to free streaming of particles). Thus any periodic variation in the particle flux de- 
tected at a distance of ~ 1 AU must arise due to in situ acceleration of particles. It is 
possible that g-modes excite Alfven waves at the solar surface and that these waves in turn 
accelerate particles locally where they are observed. However, it is unclear if the observed 
magnetic fleld variation can lead to the variation of the high energy particle flux reported 
by TML. This is a complex problem which should be looked into carefully. 

If we assume that perturbations associated with g-modes are imprinted on the solar 
wind magnetic fleld near the sun, which are carried away at the wind speed (away from 
the sun, the Alfven speed is much less than the wind speed), then the phase shift observed 
in the variation of the magnetic fleld (and thus in the particle flux variation), at a distance 
(/, over the course of a mode period, P, is ~ 2TT[d/ PVu;){SVu; /Vu;)^ where is the mean 
thermal wind speed and SV^ is the random variation of the wind speed over a time period 
P. Thus, the phase coherence time of a mode, observed in the solar wind, is estimated to 
be ~ P{PVw /(^^{Vw /SVwY ■ phase coherence time of a mode of period ~ 2.6 hours 

(frequency 105/iIIz), due to random variations in the wind speed of 1% in 2.6 hours, should 
therefore be of order a day. The coherence time observed by TML (flg. 2 of their paper) 
appears, however, to be about a year. 

3. Velocity amplitude of g-modes due to turbulent excitation 

A number of people have investigated the linear stability of solar g-modes (e.g. Dilke 
& Gough 1972; Rosenbluth & Bahcall 1973; Christensen-Dalsgaard et al. 1974; Shibahashi 
et al. 1975; Boury et al. 1975; Saio 1980). All of these investigations flnd that g-modes of 
radial-order (n) greater than 3 are stable. However, there is no general agreement about 
the stability of low order modes (n < 3). The reason for this is that there is a delicate 
balance between driving and damping for very low order modes, whereas higher order 
modes are stabilized by a rapid increase in damping with n. Thus it is almost certain that 
the high order g-modes reported by Thomson et al. can not be excited due to overstability. 

Overstability, however, is not ruled out for very low order g-modes (n < 3 & > 
150/iHz). If overstable, the g-mode amplitude will increase exponentially with time until 
nonlinear effects become important and saturate their growth. Kumar and Goodman 
(1995) have recently investigated 3-mode parametric interaction, a very efflcient nonlinear 
process. Using their results we flnd that the low order overstable g-modes in the sun will 
attain an energy of at least 10^^ erg before they are limited by nonlinearities. The velocity 
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at the solar surface corresponding to this energy is ~ 10^ cm s~"^, which is an order of 
magnitude larger than the observational limit of Palle (1991). Thus even low order g-modes 
of frequency greater than about 150 /iHz are unlikely to be overstable. 

Stable g-modes can be excited by the same process that is responsible for exciting the 
p-modes, i.e. interaction with the turbulent convection. We follow the work of Goldreich, 
Murray and Kumar (1994) to estimate the energy and surface velocity amplitude of low 
degree g-modes. The excitation is due to the fluctuating Reynolds' stress and the theory 
has been calibrated so that it flts the energy input rate in the p-modes over the entire 
observed frequency range. We also estimate the turbulent and radiative damping of g- 
modes. These damping time scales are of the same order for modes of low n (~ 10^ years), 
but the radiative damping rate increases with n as whereas the turbulent damping rate 
is a weak function of mode order. Thus radiative damping dominates at n greater than 
a few. Using the energy input rate and the damping rates we calculate the mode energy 
and the surface velocity amplitude for low degree g-modes. 

Figure 1 shows the surface velocity amplitude as a function of frequency for low degree 
solar g-modes. Note that the surface velocity amplitude is about 10~^ cm for low order 
g-modes and about 10~^ cm for moderate to high order modes of ^=1, which is a factor 
of about 10"^ smaller than the amplitude inferred from Thomson et al.'s observation. We 
also flnd that the surface velocity amplitude falls off rapidly with d making detection of 
high degree modes > 3) extremely unlikely. 

The ratio of the horizontal and radial velocity amplitudes, at the solar surface, can 
easily be shown to be equal to [d{d + l)Y GMq I {R%u'^) ^ 0.13[i{i + 1)]^/^P^ where Pi 
is the mode period in units of one hour. Thus, it is best to make Doppler measurements 
near the solar limb to search for low frequency g-modes. 

4. Conclusion 

Thomson et al. (1995) suggest that there are periodic variations in the high energy 
particle flux of the solar wind which they attribute to solar g-modes. They also suggest 
that g-modes are present in the solar magnetic fleld data. Since we expect there to be 
a variation of the thermal particle flux associated with the magnetic fleld variation, the 
oscillations should also be present in the flux of thermal particles, a testable prediction. 
Using the observed particle flux and magnetic fleld variations associated with g-modes, as 
reported by Thomson et al. (1995), we have estimated individual g-mode energies and flnd 
them to be at least 10^^ erg (for comparison, the total energy in all of the solar p-modes 
is about 10^^ erg). The corresponding surface velocity amplitude for a low order g-mode 
is at least 50 cm s~"^, which is larger than the observational upper limit (Palle, 1991, and 
references therein) by about a factor of 10. 

The phase shift experienced by a g-mode in the solar wind, as a result of a 1% variation 
in the thermal wind speed over the course of the mode period, should lead to an observed 
mode coherence time of order a day (for a 105/iHz mode), substantially shorter than that 
observed by TML (see §2). 
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High order solar g-modes (n > 3, or < 150/iHz) are found to be stable by all of 
the published linear stability calculations. These g-modes can, however, be excited by 
turbulent stresses in the convection zone. The resulting velocity amplitude at the solar 
surface is estimated to be about 10~^ cm for low order, low degree modes (our calcu- 
lations of p-mode surface velocity amplitudes due to turbulent excitation fit the observed 
velocity spectrum well, and thus the error in the predicted g-mode velocity amplitudes 
is unlikely to be greater than an order of magnitude). This is about four orders of mag- 
nitude smaller than the amplitude implied by the observations of Thomson et al. The 
ground based network of telescopes designed to observe the solar oscillations continuously 
(GONG) may be able to detect g-mode pulsation down to about 1 cm s~"^, and SOHO 
is expected to have a detection threshold of about 10""^ cm (Gabriel et. al, 1995). 
Thus g-modes, if stochastically excited, may be at the threshold of detectability of the 
forthcoming oscillation instruments. The expected long life time of g-modes (~ 10^ years 
for low order g-modes, and 10^ years for high order modes) should help in detecting these 
small amplitude pulsations. 

Acknowledgment: We thank John Belcher, Tom Duvall, Peter Goldreich, Peter Stur- 
rock, David Thomson, Louis Lanzerotti, and the referee for useful comments. 
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FIG. 1. — Magnitude of the surface velocity amplitude as a function of frequency for low 
degree solar g-modes excited by coupling with turbulent convection. The surface velocity 
amplitude falls off strongly with increasing d making detection of high d modes unlikely. 
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